Proto-Neutron Star Cooling with Convection: The Effect of the Symmetry Energy 



O 

o 

Q 



X 

o 



> 

m 
cn 
o 



X 



L.F. Roberts/ G. Shen,^ V. Cirigliano,^ J.A. Pons,^ S. Reddy,^'^ and S.E. Woosley^ 

'Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064 USA 
^Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544 
^ Departament de Ftsica Aplicada, Universitat dAlacant, Ap. Correus 99, 03080 Alacant, Spain 
^Institute for Nuclear Theory, University of Washington, Seattle, Washington 98195 

We model neutrino emission from a newly born neutron star subsequent to a supernova explosion 
to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at 
liigh baryon density. We find the time period and spatial extent over which convection operates 
is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When 
convection ends within the proto-neutron star, there is a break in the predicted neutrino emission 
that may be clearly observable. 
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The hot and dense proto-neutron star (PNS) born sub- 
sequent to core-collapse in a type II supernova explo- 
sion is an intense source of neutrinos of all flavors. It 
emits the 3 — 5 x 10^'^ ergs of gravitational binding en- 
ergy gained during collapse as neutrino radiation on a 
time scale of tens of seconds as it contracts, becomes 
increasingly neutron-rich and cools. Cooling of the PNS 
and the concomitant neutrino emission are driven by neu- 
trino diffusion and convection along the lepton number 
and entropy gradients left behind within the PNS after 
core-bounce, where the matter density and temperature 
are in the range p = 2 — 6 x lO^'^ g/cve? and T = 5 — 40 
MeV, respectively [H-Q- While the supernova explosion 
mechanism and associated fall back of material are ex- 
pected to influence the neutrino emission at early time 
(i.e. t < 1 s post bounce) the late time neutrino signal is 
shaped by the properties of the PNS, such as the nuclear 
equation of state (EoS) , neutrino opacities in dense mat- 
ter, and other microphysical properties that affect the 
cooling timescaie by influencing either neutrino diffusion 
or convection [3, Ull] . 

Here, we present one-dimensional hydrostatic models 
of PNS evolution out to late times for two EoSs. We 
include approximate convective transport along with dif- 
fusive neutrino transport, both consistent with the un- 
derlying EoS. This allows us to gauge the importance of 
convection and the effect of medium modifications to the 
neutrino opacities in dense matter to the temporal char- 
acteristics of the neutrino signal. The basic framework 
for PNS evolution is similar to that described in ex- 
cept that a treatment of convection is included. We find 
that the behavior of the nuclear symmetry energy (which 
is a measure of the energy difference between dense neu- 
tron matter and symmetric nuclear matter) at high den- 
sity significantly influences convection and thereby affects 
the observable neutrino signal. 

Large scale convective overturn of material will directly 
transport energy and lepton number in the PNS and al- 
ters the gradients along which neutrinos diffuse, thereby 
strongly affecting the neutrino signal accompanying PNS 



formation. It has long been recognized that the outer 
PNS mantle is unstable to convection soon after the pas- 
sage of the supernova shock, due to negative entropy gra- 
dients [13]. This early period of instability beneath the 
neutrino spheres has been studied extensively in both one 
and two dimensions, with the hope that it could increase 
the neutrino luminosities enough to lead to a successful 
explosion [2, 11-15|. Although the role of convection at 
late times was studied in Refs. this work is the first 

attempt at exploring its connection to the underlying mi- 
crophysics and the interplay between modified opacities 
and convection in shaping the observable neutrino signal. 

We model convection as a diffusive process described 
by standard time dependent mixing length theory 0] ■ In 
Ref. [isl it was shown that hydrodynamic simulations in 
two-dimensions were well reproduced by a simple mixing 
scheme in one-dimension, suggesting that this is a rea- 
sonable approximation. Further, since convection pro- 
ceeds efficiently throughout our simulations, our results 
are not particularly sensitive to the chosen mixing length. 
The linear growth rates are obtained from the standard 
Ledoux stability analysis as 

Lo^ = -^{1sVHs)+1y,VHYl)), (1) 
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and g, P, s, ns, and Yl are the local acceleration due to 
gravity, the pressure, entropy per baryon, baryon num- 
ber density and the fraction of leptons in dense matter, 
respectively. Convective instability sets in for > 0. 
The form of the growth rate clearly implies that there 
will be a strong interplay between the nuclear EoS and 
the patterns of convection within the PNS. 

In addition, the PNS may be subject to doubly diffu- 
sive instabilities due to the lateral transport of compo- 
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FIG. 1: The symmetry energy as function of density for the 
lU-FSU and GM3 EoSs. Inset: noiSsym versus Esym at nuclear 
saturation density, for lU-FSU (circle), GM3 (square), and 
QMC (diamonds). The shaded regions correspond to various 
experimental constraints taken from Ref. [l^ . 



instabilities and convection enhanced the neutrino lumi- 
nosity to successfully power a neutrino driven explosion. 
However, more recent two dimensional studies found no 



evidence of these doubly diffusive instabilities [IJ, Il5l |. 
Because of this and the increased complexity of treating 
the doubly diffusive instabilities, we do not include them 
in our study. 

The EoS and neutrino interaction rates are modeled 
using a relativistic mean field (RMF) model of nuclear 
interactions. We adopt a non-linear generalization of 
the original Walecka model described in [l^. Here, 
the nucleon-nucleon interaction energy is calculated in 
the mean field approximation using effective interactions, 
which are tuned to reproduce gross observed properties 
of nuclei and empirical properties of symmetric nuclear 
matter at saturation density. Although these empiri- 
cal constraints provide valuable guidance to constrain 
aspects of the symmetric nuclear EoS at nuclear den- 
sities, the experimental constraints on the properties of 
neutron-rich matter are relatively weak. The difference 
between the energy of symmetric matter (equal num- 
bers of neutrons and protons) and pure neutron matter is 
called the symmetry energy, -Esym('^B), and is defined by 

E{nB,Xp) = E{nB,Xp = 1/2) + Esym{nB)S^ H . Here, 

d — {I — 2xp) and E{nB, Xp) is the energy per particle of 
uniform matter composed of neutrons and protons with 
total baryon density ub and proton fraction Xp . In charge 
neutral matter Xp — where is the electron fraction. 
Various experimental probes of the nuclear symmetry en- 
ergy and its density dependence in nuclei and heavy-ion 
collisions are actively being pursued in terrestrial exper- 
iments, but are yet to yield strong constraints. These 
constraints are shown in the inset in Fig. [1] and are dis- 
cussed in Refs. Quantum Monte Carlo (QMC) 
results are also shown in the inset in Fig. [TJ The linear 
correlation between i?sym and E'^^^-^ in the QMC results 



is obtained by varying values of the poorly known three- 
neutron interaction [20|. 

Recent work has shown that the derivative of the sym- 
metry energy with respect to density, denoted as E'^^j^ = 
dEsyml dnB, plays a crucial role both in the terrestrial 
context where it affects the neutron density distribution 
in neutron-rich nuclei and in astrophysics where it affects 
the structure and thermal evolution of neutron stars (for 
a recent review see Ref. [131 )• The pressure of neutron 
matter at saturation density, Pneutron(fiQ) = n^Ei^,^, in- 
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fluences the radii of cold neutron stars 
rich nuclei, the neutron-skin thickness is also sensitive 
to £'gyjjj(/9o), so that there exists a linear correlation be- 
tween the neutron-skin thickness and neutron star radius 
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To study the sensitivity of PNS evolution to the nu- 
clear symmetry energy we employ two RMF models with 
different predictions for EL^{po). The first EoS is the 
lU-FSU EoS taken from [19| , which includes a non-linear 
coupling between the vector and iso-vector mesons that 
allows the symmetry energy to be tuned at high den- 
sity. The second EoS employed is the GM3 parameter 
set, where non-linear coupling of the vector meson fields 
is neglected [2l|. The symmetry energy as a function of 
density is shown in Fig. [T]for the two EoS. The inset in 
Fig.[T]shows current theoretical estimates and experimen- 
tal constraints on E'sym and hqE'^^^ at nuclear density. 

In the rest of this letter, we demonstrate that E!,y^(pQ) 
plays a role in stabilizing PNS convection at late times 
and thereby directly affects the PNS neutrino signal. The 
logarithmic derivatives 7s and are always positive, so 
that negative entropy gradients always provide a destabi- 
lizing influence. For given entropy and lepton gradients, 
stability is then determined by the ratio JYl/Is- The 
sign and magnitude of is strongly influenced by the 
density dependence of the nuclear asymmetry energy, so 
that negative gradients in lepton number can be either 
stabilizing or destabilizing and the degree to which they 
are stabilizing varies from EoS to EoS. To clarify this we 
note that at T = and when the neutrino contribution 
to the pressure is small 
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which is a reasonable approximation to the finite tem- 
perature result. The first term comes from the electron 
contribution to the pressure, while the second term is 
due to nucleons and is negative since both the Fermi and 
interaction energies favor a symmetric state. For high 
densities and low electron fractions, for realistic values 
of £'syni, this leads to negative 7y^. Therefore, a larger 
iJgyjjj leads to negative lepton gradients in the PNS pro- 
viding a larger stabilizing influence. E'^y^ also partially 
determines the equilibrium value of Yg, which can alter 
the value of 7^^^ , but this is a smaller effect. In our nu- 
merical PNS simulations this effect is accounted for. In 
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FIG. 2: Evolution of the entropy (top panels) and lepton frac- 
tion (bottom panels) in a 1.6Mq rest mass PNS for the GM3 
EoS (left panels) and the lU-FSU EoS (right panels). The 
grayed regions are convectively unstable. The labels corre- 
spond to the model times in seconds. 



contrast, the properties of zero temperature nuclear mat- 
ter should have significantly less effect on the behavior 
of 7s. We find the variation of 7s to be significantly less 
than that of jyl in ttie two EoSs considered here, which 
is consistent with expectations. 

Medium modifications of the neutrino interaction rates 
also infiuence PNS evolution. These effects have been 
investigated in earlier work where neutrino scattering 
and absorption rates on nucleons and leptons were calcu- 
lated within the relativistic random phase approximation 
(RPA) [2^. Here, the effects due to strong and electro- 
magnetic correlations, degeneracy, and relativistic cur- 
rents are included in a similar way. Further, the effective 
interactions we employ in the RPA are obtained from the 
underlying model for the EoS to ensure that the corre- 
lation functions which determine the neutrino scattering 
kernels satisfy the generalized compressibility sum-rules 
consistent with the EoS. Because the axial portion of 
the response accounts for the majority of the neutrino 
scattering cross-section, the presence of an axial inter- 
action can strongly affect the total neutrino mean free 
path. To take this into account, we introduce an effec- 
tive short-range interaction in the spin channel through 
the Migdal parameter, g' [i^. The strength of this in- 
teraction is tuned to reproduce the spin-susceptibility of 
neutron matter obtained from microscopic calculations 
27| . For densities above nuclear saturation, the RPA 
causes a significant enhancement of the mean free path 
relative to the mean field approximation due to the re- 
pulsive nature of the nuclear interaction at high density. 

We now consider the evolution of the internal structure 
of the PNS with convection and varying prescriptions for 
the opacities. In Fig. [21 the evolutions of the entropy 
and lepton fraction for the two equations of state are 
shown. Over the first second in both models, convec- 
tion smoothes the entropy and lepton gradients in the 



outer regions to a state close to neutral buoyancy. GM3 
has a slightly steeper entropy gradient, in part due to 
its larger E'^y^, than lU-FSU. This partially accounts 
for the slightly larger neutrino luminosity at early times 
for GM3. As time progresses, convection steadily digs 
deeper into the core of the PNS. For both EoSs, con- 
vection proceeds all the way to the core by 15 seconds 
into the simulation, but it lasts in the interior regions 
for a much longer period of time for lU-FSU resulting 
in more rapid lepton depletion in the core. The exact 
details of how convection proceeds depend on the initial 
conditions of the PNS and the behavior of 7s and jyl 
for a given EoS. For the conditions encountered in the 
PNS, the variation in 'JYl between the two EoSs we em- 
ploy is significantly larger than the variation of 7s . More 
important to the neutrino signal accompanying PNS for- 
mation, in GM3 convection ceases in the mantle by ^ 5 
seconds, whereas convection in the mantle proceeds until 
~ 12 seconds in lU-FSU. This difference is mainly driven 
by the difference in E!.y^ between the two EoSs. As the 
mantle contracts, the second term in Eq.([2]) becomes in- 
creasingly dominant and is eventually able to stabilize 
convection. Qualitatively, increasing E!.y^ will shut-off 
convection at an earlier time. 

The depth to which convection penetrates in the core 
and how long convection proceeds in the core are depen- 
dent upon the opacities as well as the EoS. When only 
mean field effects on the opacities are considered (i.e. 
when the neutrino mean free path is shorter), convection 
does not proceed all the way to the center of the PNS in 
the GM3 models. When RPA effects are included, con- 
vection does proceed to the central regions of the core. 
An increased diffusion rate allows the core to heat up and 
deleptonize more rapidly, thereby decreasing the stabi- 
lizing lepton gradients and increasing the de-stabilizing 
entropy gradients. 

Of course, variations in the convective evolution of the 
PNS are only interesting to the extent they are poten- 
tially observable in the neutrino emission from a nearby 
supernova. In Fig. |3] the expected neutrino count rates 
for a detector similar to Super Kamiokande-III are shown 
for a number of PNS cooling models. We have assumed 
a threshold energy of 7.5 MeV, a detector mass of 50 kt, 
a detector efficiency above threshold of unity [2S] , and a 
distance of 10 kpc to the supernova. Equipartition has 
been assumed between neutrino flavors and the spectral 
temperatures are calculated from the average radius of 
neutrino decoupling in the PNS. The neutrino luminosi- 
ties and average energies for two models are included as 
supplemental data. 

Both the GM3 and lU-FSU EoSs show enhanced lumi- 
nosities at early times relative to the models not includ- 
ing convection. There is only a small difference between 
the two equations of state at low (sub-nuclear) density, 
so differences prior to one second are small. The neu- 
trino count rate is increased by about 30% relative to the 
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FIG. 3: Count rates as a function of time for a number of 
1.6Mq PNS models with and without convection. The black 
line is for neutrino opacities calculated in the mean field ap- 
proximation, while all the other lines are for models that use 
RPA opacities with g' = 0.6. The inset plot shows the inte- 
grated number of counts from 0.1 s to 1 s divided by the total 
number of counts for t > 0.1 second on the horizontal axis, 
and the number of counts for t > 3 seconds on divided by the 
total number of counts for t > 0.1 second. The stars corre- 
spond to the lU-FSU EoS and the circles to the GM3 EoS. 
Symbol sizes correspond to various neutron star rest masses 
ranging from 1.2Mq to 2.1Mq. Colors correspond to different 
values of the Migdal parameter, g' . 



models that do not include convection. This is reason- 
ably consistent with the e arly time enhancement seen in 
multi-dimensional models [15| . After a second, the count 
rates between the two EoSs begin to diverge. The most 
obvious feature in the count rate for GM3 appears at ^ 3 
seconds, which is coincident with the end of convection 
in the mantle. For the lU-FSU EoS, the break is also 
at the time at which mantle convection ends (~ 10 sec- 
onds), although it is hard to distinguish from the point 
at which the PNS becomes optically thin. As was ar- 
gued previously, the position of this break reflects the 
density dependence of the nuclear symmetry energy at 
riB > no and therefore provides a direct observable of 
the properties of nuclear matter in the PNS neutrino sig- 
nal. Although core convection does not seem to affect the 
break, it may impact the subsequent cooling timescale. 

In the inset in Fig. [3l integrated neutrino counts over 
two time windows are shown for a number of PNS masses. 
There is a clear separation between the two EoSs inde- 
pendent of mass. The time of the convective break cre- 
ates this separation. This illustrates that this diagnostic 
of the symmetry energy does not require an accurate de- 
termination of the PNS mass. 

The inclusion of nucleon correlations through the RPA 
begins to significantly affect the neutrino emission after 
about three seconds. Initially, the luminosities are in- 
creased as energy and lepton number are able to more 
rapidly diffuse out of the core, but at later times the 



neutrino signal is significantly reduced and drops below 
the detectable threshold at an earlier time. 

In summary, using a self-consistent model for the PNS 
core physics, we find that the late time neutrino signal 
from a core collapse supernova is likely to contain a di- 
rect diagnostic of the nuclear symmetry energy at high 
density. With current neutrino detectors, these effects 
should be readily discernible in the neutrino light curve 
of a single nearby supernova. 
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